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Effect of pair-breaking on mesoscopic persistent currents well above the 
superconducting transition temperature 
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We consider the mesoscopic normal persistent current (PC) in a very low-temperature supercon- 
ductor with a bare transition temperature T° much smaller than the Thouless energy E c . We show 
that in a rather broad range of pair-breaking strength, T,? < h/r s < E c , the transition temperature 
is renormalized to zero, but the PC is hardly affected. This may provide an explanation for the 
magnitude of the average PC's in the noble metals, as well as a way to determine their T°'s. 

PACS numbers: 74.78.Na, 73.23.Ra, 74.40.+k, 74.25.Ha 



Introduction. The magnitude of the equilibrium aver- 
aged persistent currents (PC's) [l|, [2J in normal metals 
has been a long-standing puzzle. Experiments [3, U, M, 
produce a current larger by at least two orders of mag- 
nitude than the theoretical prediction for nonintcracting 
electrons jy, @, Q and seem to indicate that the low-flux 
response is diamagnctic. The average PC of a diffusive 
system with interactions was calculated first in this con- 
nection [9j] in Refs. [lfj and [ll|. The Resulting PC was 



1 --B- 

0.5 





l/l(s=0) 




0.1 



10 



found to be much larger than that of a noninteracting 
system, but nevertheless not large enough to explain the 
experiments. 

Repulsive electron-electron interactions [ll( result in a 
paramagnetic response (at small magnetic fluxes) whose 
magnitude is smaller than the experiment by about a 
factor of five. This disagreement is due to the down- 
ward renormalization of the interaction [12L Il3| . Attrac- 
tive interactions [10( result in a diamagnetic response, 
whose magnitude (due to the very low superconducting 
transition temperature), is again smaller by a factor of 
order five than the measured one. This is in spite of 
the renormalization upward of the attractive interaction. 
Attractive interactions, at low energies, imply (with no 
pair-breaking) a transition into a superconducting state, 
and the PC of such an interacting system depends on its 
transition temperature. These temperatures are very low 
14 1 for the noble metals used in the PC experiments - 
hence the too small predicted values for the PC. 

Here we consider attractive interactions. We show that 
the presence of a very small amount of pair breakers, 
e.g., magnetic impurities (which seem to be very diffi- 
cult to avoid in these metals |la|), may change the pic- 
ture profoundly. Obviously one may consider other pair- 
breakers, such as a two-level systems [16| or simply a 
magnetic field [l7(. In this letter we treat specifically 
the case of magnetic impurities. We find that within a 
significant range of the pair-breaking strength, the mag- 
netic impurities suppress the transition temperature down 
to immeasurable values, leaving concomitantly the PC al- 



FIG. 1: The first flux harmonic [m — 1, see Eq. [2], in 
units of I(s = 0), of the PC at T = E c (full line) and 
T c /T c (dashed line) as functions of the pair-breaking strength, 
s — 1/(ttT°Ts), displayed on a logarithmic scale. 



most unchanged. The physical reason for this remarkable 
observation is that the PC is determined by the interac- 
tion on the scale of the Thouless energy E c = HD/L 2 
(~ 20mK for a typical experimental system), while the 
bare transition temperature, T c °, is much smaller. (The 
circumference of the ring is denoted by L and D is the 
diffusion coefficient.) This gives rise to a rather wide 
range of pair-breaking strengths, presented here by the 
spin-scattering time r s , 



r c ° < h/ Ts < E c , 



(1) 



in which the actual transition temperature T c will drop 
to zero [18| . but the PC will be hardly affected. As a 
result, it is the bare transition temperature of the system 
without the magnetic impurities, Xj?, as opposed to T c , 
which dominates the expression for PC, see Fig. Q] We 
concentrate here on the experimental results of Ref. y 
[191 ]. In order to explain them, it is necessary to assume 
a T~ in the lmK range for copper. Our basic assertion is 
that this may indeed be the correct order of magnitude of 
T° for ideally clean copper, but that it is knocked down to 
zero or to a very low value by a minute, < ppm, amount 
of unwanted ]laJ pair-breakers. We emphasize, however, 
that our result concerning the fundamentally different 
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FIG. 2: The first flux harmonic of the PC in units of 
I* = —eE c as a function of the temperature, for two val- 
ues of T°/E c and several values of s. Keeping, at T < E c , up 
to the 100 lowest values of \u\, was necessary for convergence. 
Note that the s = curve in the upper panel is valid only for 
T/T c > 1 + Gi, where Gi is the Ginzburg parameter (Gi ~ 0.1 
for the samples of Ref. (3). 



sensitivities of T c and PC to pair-breaking in the range 
given by Eq. ([1]), remains valid regardless of the situation 
in specific materials. The Kondo screening of the spins is 
not considered here. Other effects of magnetic impurities 
have previously been considered in Ref. |20j. 
Results. The expression we obtained for the PC in a 
diffusive ring with magnetic impurities can be expressed 
as a sum over the harmonics of the magnetic flux through 
the ring t/>, in units of the flux quantum h/e, 
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1). Here v denotes the bosonic Matsubara 



(using h 

frequency 2l|, 'F and \&' are the digamma function and 
its derivative, and T is the temperature. Our expression 
^ generalizes the result of Ref. [lC] for the case where 
spin-scattering is present: the Matsubara frequency \v\ 
is shifted by 2/r s . However, the superconducting tran- 
sition temperature (which appears formally in the de- 
nominator of the integrand) is not the one modified by 
the pair-breakers, but retains its bare (magnetic impuri- 
ties free) value. Interestingly enough, it follows that by 
measuring the PC one may determine T® (which would 
be directly measurable only if all low-temperature pair 
breaking could be eliminated). 

In Fig. [2] the PC is plotted using Eq. fl2|). At the 
critical pair-breaking time l/r s ~ T®, corresponding to 



s = 1/itt s T® ~ 1/7T, the transition temperature vanishes 
[18(|, while the PC is hardly affected. The measured PC 
in the copper samples of Ref. y is I(T < E c ) ~ —eE c . 
The curve with s = 1 in the upper panel, taken with 
T c ° = l.5mK and E c = \hmK (the value for the samples 
of Ref. y) gives a PC lower by only 25%. A better fit is 
possible by changing the parameters somewhat, but we 
do not regard this as crucial at the present stage. Like- 
wise, we can qualitatively explain the result of Ref. |4J. 
The high frequency results of Ref. |5j require a sepa- 
rate discussion [5j. The PC is reduced significantly once 
1/t s > E c , or L s = ^/I>F S < L. For T?/E c = 0.1 (0.01), 
the condition for E c t s ~ 1 is s — 10 (100). 
Derivation. For completeness, we outline below the 
derivation of the PC in the presence of magnetic scatter- 
ing [lj|. The PC, Eq. @, is obtained by differentiating 
the free energy with respect to the flux. Our system is 
described by the Hamiltonian [18[ 

H = fdr(^i(r) [(Ho + «i(r))<5 Q7 + u 2 (r)S ■ <r^] V 7 (r) 
- §V£(r)^(rW 7 (r)tf a (r)) , (3) 

in which the last term is the attractive interaction, of cou- 
pling g. The spin components are a and 7, a is the vector 
of the Pauli matrices, and Ho = (— i V — eA) 2 /2m — /i (/i 
is the chemical potential and A is the vector potential de- 
scribing the flux through the ring). The scattering, both 
nonmagnetic and magnetic, is assumed to result from N{ 
point-like impurities, such that 



ui(r) + u 2 (r)S • er 



N ' 1 

J^(s(r - Ri) - -) (ui + u 2 S Rt ■ a) 



(4) 



where V is the system volume. In averaging over the im- 
purity disorder one assumes that the impurity locations, 
Ri, are random, and so are their classical spins, such that 



'R, 



0, and <S R • Sj, > = SijS(S + l). 



The partition function, Z, is calculated by the method 
of Feynman path integrals [22J , combined with the Grass- 
man algebra of many-body fcrmionic coherent states in 
terms of the variables ip a (r,r) (^(r, r)). Introducing 
the bosonic fields A(r,r) via the Hubbard-Stratonovich 
transformation leads to the partition function Z = 



J D(i/>(t, r), </>(r, r))D(A(r, r), A* (r, T ))t 



with 



dr f e dr (\^Dl 

Jo v 9 
i§(r,r)G r -i ;T , T *(r,r)) , 



(5) 



where \& = (ipijipi^tpi, tpi). The inverse Green function 
G _1 (at equal positions r and equal imaginary times r) 
is given by 



G~l _ , = 
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where h„^ = Hq(±A) + u\ + sgn(a)S z ii2, and S± = 
(S x ±iS y )/2. 

The integration over the fermionic part of the action 
© yields 



Z = JD(A(r,T),A*(r, T )) 



/I 



exp(-Trln(/3G" 



- / dv 



(It 



A(iy 



(7) 



In order to treat the boson fields A, we expand 
Trln(/3G _1 ) up to second order in A. This expansion is 
valid for temperatures well above the transition temper- 
ature, and, strictly speaking, above the Ginzburg critical 
region. The zeroth order is omitted as it leads to the 
tiny magnitude PC of noninteracting, grand-canonical, 
normal metal rings [6| . The result in Fourier space reads 
(the dependence on the magnetic flux is specified below) 



TrlnOSG- 1 ) 



2 1 1,1 



J2 J2 Tr G , p(k 1 + qi,k 2 + q2,w + ^)A(q 2 ,z/)G , h(k 2 ,k 1 ,-tj)A T (q 1 ,; 
qi,q2,^ ki,k 2 ,w 



(8) 



r 



The resulting expression for the partition function may 
be simplified considerably. Firstly, the terms that survive 
the disorder-average in the sum of Eq. ([8]) are those for 
which [23j qi = q 2 . Secondly, the particle and the hole 
Green functions, G p and Gh, [see Eq. ©] are related, 



G h (k,k';^) = -G P (-k,-k» , 



(9) 



where the superscript t denotes the transposed. Carrying 
out the integration in Eq. (J7J), 






(10) 



where 7V(0) denotes the extensive density of states at the 
Fermi level. The polarization is 



n(q, v) = - Y^ £ oryKuory(<l, v) 



(11) 



with 



K uaj (q, v) = 22 ( G aa> (ki + q, k 2 + q, lu + v) 

ki,k2 

x e Q / 7 ,G 7 y(-k 2 ,-ki,-w)) . (12) 

Here e is the anti-symmetric tensor, e aa = 0, and e^ = 
— £j.t = 1, and G denotes the particle Green function. 

In Ref. lla K(0, 0) was calculated using a Dyson equa- 
tion. We generalize their calculation to obtain K(q, v) 



and consequently the polarization becomes [19( 
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Here cod is the cutoff frequency on the attractive inter- 
action, and the pair-breaking time r s is given by 



1 

T. 



2irN(0)NiS(S+l)ul 



(14) 



The transition temperature of the system in the ab- 
sence of pair breakers, T°, is obtained from the q — 0, v — 
pole of Z, upon setting 1/t s = 0, 



— =*ri 

9 U(0) V2 



up 
2ttT? 



< 



(15) 



(Note that the same procedure in the 'presence of the 
pair breaking reproduces the decrease in the transition 
temperature T c , as found in Ref. UM-) Since u>d ^> T^},T 
we may use the asymptotic expansion of the digamma 
function. In this way we obtain 



n[-Q 



q.// 

/! Dq 2 + |H+2/r 8 
V2 4ttT 
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(16) 



Finally, the PC is given by I = (e/h) dT\nZ/dcj). In 
our ring geometry, the flux enters the longitudinal com- 
ponent, qu, of the vector q as 



q\\ = -^{n + 2(j)) , 



(17) 



where n is an integer. Only the zero transverse momen- 
tum contributes significantly to the current. Our result 
(|2|) is obtained upon inserting Eq. (fl6|) into the defini- 
tion of the current and employing the Poisson summa- 
tion formula. It then follows from Eq. @ that values of 
t s which are detrimental to T c , may hardly affect the PC 
(see Fig. H}. 

We conclude by further explaining the physical argu- 
ment behind our result. Very roughly, the renormaliza- 
tion of the dimcnsionlcss attractive interaction A (> 0) 
from a higher frequency scale w> to a lower one, w<, is 
given by A(w<) = 1 _ A( ^ ) f n ) ( ^ ) • At T c ° and l/r s = 0, 

the attractive interaction should diverge. Using this 
to eliminate \(ujd) (= gN(0)/V), we obtain that for 
T° < w « U) D , A(w) - l/7n(w/T c °), which around the 
Thouless scale, is close to the value found in Ref. lid . 
The pair-breaking stops the renormalization at l/r s , but 
does not significantly change the interaction on the much 
larger scale of E c . Our prediction can also be tested with 
very small rings made of known low T c superconductors. 
We point out that the mechanism suggested by 
Kravtsov and Altshuler 24j , relating extrinsic dephasing 
to an enhanced PC, is different than ours, since it relies 
on the rectification of the noise. 
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